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Abstract: We consider tunneling between 2 symmetric potential wells for 
a 2-d Schrodinger operator, in the case of eigenvalues associated with with 
quasi=modes supported on KAM or Birkhoff tori. 

0. Introduction. We consider here tunneling between 2 symmetric potential 

wells for a 2-d Schrodinger operator P — —h^A + V in the limit h ^ 0, near 
some energy level Eq = close to the non degenerate minima of V. 

Tunneling is a difficult problem that has exercised so far many subtle and 
ingenious strategies ; at least, computing tunneling rates involves various sce- 
narios which depend on the details of the dynamics, ranging from integrable 
or quasi-integrable systems, to ergodic or chaotic ones (see [W], and [Cr] for a 
recent review. ) 

As a general rule, the energy shift (or splitting of eigenvalues) is related 
to the so called Agmon distance So{E) between the wells, associated with the 
degenerate, conformal metric ds'^ = (^ ~ E)+dx^ that measures the life-span 
(instanton) of the particle in the classically forbidden region V{x) > E. 

Much is known in the 1-d case, even for excited states, or in several di- 
mensions for the lowest eigenvalues. For general wells, there is the following 
equivalence [Ma] : the corresponding normalized eigenfunctions are non expo- 
nentially small (i.e. for all e > 0, larger, in local norm, than a constant 
times e"^/'', < h < /ig, ) where minimal geodesies, connecting the 2 wells, 
meet their boundary, if and only if the splitting is non exponentially small with 
respect to e'^"^^'^/'^. 

We study the special case of splitting of eigenvalues associated with quasi- 
modes supported on KAM or Birkhoff tori ; our goal is to compute tunneling 
rates for a large family of such eigenvalues, which we shall call a spectral tunnel 
series. 

1. Tori and quasi- modes. Let us consider for a moment the case of a single 

well around Uq = (so that we can ignore interaction with the other well, ) and 
let Po{x,S,) = -|- Xfxf + X2X2 be the quadratic part of the (smooth) classical 
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hamiltonian p{x,£^) near (quadratic approximation). If the frequencies Aj's 
are rationnally independent, then Birkhoff's theorem tells us that the orbits of 
p near the fixed point (x, ^) = with energy E are quasi-periodic, in the sense 
that they arc confined within quasi-invariant tori (the Birkhoff tori,) over a 
time scale 0{E~°°). Whenever the system is non integrable, most of these tori, 
will be destroyed and replaced by chaotic regions ; however (under a suitable 
a diophantinc condition on the frequencies, ) the so-called KAM tori, whose 
collection form a Cantor set, eventually survive. 

Such lagrangian (or possibly isotropic) integral manifolds Af support quasi- 
modes, whose energies arc given by the Einstein-Brillouin-Kellcr (EBK) quanti- 
zation rule : if H{l) = 2Ai(-i -|- 2X2L2 + ■ ■ ■ is the hamiltonian expressed (asymp- 
totically) in terms of action-angle variables (t, if), then the energies of the eigen- 
states are given to first order in /i by : 

(1) E„{h) = H{La), t^ = {a + ^iy)h, a e N"^ , \a\h < Eq 

The vector 1/ = (2, • • • , 2) is Maslov index, counting the number of caustics met 
along an orbit, passing from sheet to sheet on the torus. 

For Birkhoff tori associated with energies Ch < E < , where < 5 < 1, 
EBK formula (1) can be corrected to all orders in h with an accuracy 0(/i°°) [Sj]. 
In case of KAM tori, and larger energies (in an interval independent of h, ) these 
expansions hold modulo 0{e~^^^'^^^°), for some s > 1 related to the diophantine 
condition on the ratio A1/A2 [Po]. Of course in that case, analyticity properties 
of the potential are required as in usual KAM procedures, which makes the 
analysis more subtle. But the geometry is the same. 

We have seen that tunneling rates hinge upon decay of quasi-modes at the 
edge of the well OUe = {V{x) — E}. They also depend on the decay near the 
caustics of Af , «, = <,„• 

The caustics can be viewed as a rectangle shaped fold line delimiting the 
zone of pure oscillations of the quasi-modes, and touching the boundary of the 
wells at 4 vertices, the hyperbolic umbilic points (HU) points, section of the 
torus by the plane ^ = in R**. According to (1), the edges of this rectangle 
have a size Lq. — (a + ^i^)h. 

All tori Af continue analytically in the ^ variables. Over the classically 
forbidden region, analytic continuation amounts to parametrize the orbits with 
imaginary time. It is convenient to view Af' as a multidimensional Riomann 
sheet structure, with a number of sheets corresponding to the choice of the sign 
of momentum, gluing along the caustics, and all intersecting at the HU's. Af 
is parametrized by a phase function Fy{x), where y G OUe denotes the umbilic 
and can be identified with t. This phase is complex (refiecting the oscillations of 
the corresponding quasi- mode) on all but one sheet of Af , denoted by Af , which 
lies over the classically forbidden region, and corresponds to a pure exponential 
decay of the quasi-mode (we exclude the sheet which gives exponential growth.) 

In case of Birkhoff tori, we gave a complete asymptotic expansion of the 
quasi-mode in some region of the decaying zone [KaRo] , close to, but at a finite 
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distance from the caustics. Similar expansions can be obtained closer to the 
caustics, in term of special functions of Airy type. Yet another expansion could 
be found in case of KAM tori. 

2. Continuation of action integrals in the classically forbidden re- 
gion. Another central geometric figure of the problem is the integral manifold 

of q{x,^) = —p{x,i^) passing above OUe, i-e. 

Af = {exptHg{p) : p € OUe x 0, q{p) = -E, t e R} 

This is (locally) a smooth real lagrangian submanifold, of the form ^ = WdE{x), 
X ^ Ue- where ^^(a;) = dE{x,dUE) is Agmon distance from x to BUe- Ac- 
tually, has the fibre bundle structure Ag = Uj/eac/jsT!/- Here 7^ is the 
bicharacteristic of q{x,^) at energy —E issued from BUe at the umbilic y, and 

7y =Af nAf. 

Introducing appropriate coordinate charts of hyperbolic action-angle vari- 
ables ((-', If) I— > {l, ifi') given by Birkhoff transformations, we can also view Fy{x) 
as the action ^dx computed along some path keeping constant, and varying 
tf. • Our first task is to continue x > F^(x) from Af , keeping y and E fixed. 

First we introduce some scaling factors. One difficulty throughout consists 
in the range of different scales. So let /x = VE be the characteristic size of the 
(euclidean) diameter of Ue- 

Let y = (2/1,1/2) be an umbilic, and assume that the torus Af is not "too flat" 
in a certain sense, or equivalently, that the rectangle shaped caustics is "not too 
far from being a square" . For x also close enough from jy, so that Fy [x) is still 
real, we denote by y(x) the unique point of OUe such that x G ly{x)- We have : 
Proposition 1: Fy(x) equals dsix) precisely along the geodesic 7^,. Moreover 

dE{x) - F^{x) ~ ^-Fjf (a;)|xe7,;(2/i - Vii^))^ ~ -K{x,y){yi - yi(a;))^ 

where K{x,y) > Kq > whenever disi{x, BUe) < Cfi^^^. (here dist stands for 
the euclidean distance.) Denote by Ty{x) the orthogonal projection of x onto jy. 

Using a variant of Gauss Lemma (the geodesic flow is locally a radial isometry), 
Proposition 1 shows that there are smooth level surfaces Nn{s) = {dE{x) = sp,} 
, si < s < S2, dist(^Ni^{s),dUE) ~ Z^^^^, such that 

(2) dE{x)-F^{x)r^--{x-Ty{x))\ xeN^s) 

Using that eikonal equation is satisfied by both dsix) and F^(x), estimate (2) 
continues in the large, all along jy, so far as 7y does not reach any caustics. 
This holds in particular, if 7 = T^; is a minimal d^-geodesic between -/V;u(s) 
and a fixed Xq, somewhere in between the 2 wells. • Our second task is to 

compute action from continuation of energy surfaces, i.e. by varying E (and 
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i' accordingly) but keeping x fixed. Let [x,S^{x)') = (.t, V,,.d£(.T))) G Ag, and 
{z' ,C,') = k(x,^{x)), where k is a suitable canonical transform related to the 
mapping {l',ip) i-^ {l,^'), that preserves the boundary of the well : K{y,0) = 

(y',0). 

In fact, we shall compute dsix) in a (/U-independent) neighborhood w of dU e 
from Agmon distance do{x) at energy 0, which is known to be a C°° function of 

z'- 

X. We have : Proposition 2: For x Geo, dsix) = do{x) + l', log ^ +C'(/U^). 

Typically, for a; € w, (.'• log ^ is comparable to /U^ log/x, so Proposition 2 gives 

the singularity of dsix). 

We look next how does dE{x,du!) depends on E in the large. For fixed xo 
away from oj, let Te as above, be a minimal d^-geodesic between TV = du and 
xq, parametrized with arc-lenght. We can arrange so that A'' is a level surface 
for dE- Using variations of geodesies as in [HeSj], we prove : Proposition 3: 

For all e > 0, there is a (/U- independent) neighborhood flE of T£;([0, 1 — £]), a 
/i^-neighborhood Ie of E, such that : (i) (x, E') ^ dE'{x, N) € C°°{VtE x Ie). 
(ii) ^Ie is starshaped, in the following sense : V(.t, E') G 0^; x Ie, 3! rf^j'-minimal 
geodesic joining N to x that stays in Qe- 3. The tunnel cycle. We label 

objets belonging the to left (resp. right) well with subscript L (resp. R). Extend 
Al = (Af)^ along the bicharacteristic flow of q. 

For a general, non integrable system, there is no reason for this extension 
coincides with Ar = (Af)^. However, wc say that the pair (pi, p^) G A^ x Ar, 
are in correspondance along a bicharacteristic 7 if A^ (or equivalently, because 
of symmetry. An) supports a quasi-mode, and {pl,Pr) S 7 x 7. We call the 
bicharacteristic 7 a tunnel cycle if there is a pair {pltPr) in correspondance 
along 7, with {pl,Pr) G i^d) l ^ i^d) r- Then and pn are necessarily 
umbilics, and 7 a geodesic between Ul{E) and Ur{E). A tunnel cycle will be 
called minimal \i the geodesic 7 is minimal, hence of lenght dE{UL{E), Ur{E)). 
Generically, bicharacteristics connecting pairs in correspondance and tunnel cy- 
cles are discrete sets. Moreover, pairs in correspondance, in case of Birkhoff tori, 
are only defined modulo 0{E°°), since this is the case for Af and Af. Tunnel 
cycles are exceptional, but as we shall see, there are many pairs in correspon- 
dance (belonging to different bicharacteristics) close to the umbilics. See [Gr] 
and [DoSh] for related notions. 

a) The case of a minimal tunnel cycle. The picture is the following : 
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Let ul{x, E, h) and ue.{x, E, h) be the quasi- modes associated with the um- 
biUcs PL and pn, continued beyond the symmetry axis S separating the 2 wells, 
in a neighborhood of the minimal geodesic T^;. Using Agmon estimates as in 
[HeSj], we can show that they approximate suitably the true eigenfunctions (pro- 
vided a gap condition. ) Assume for simplicity there is just one such minimal 
geodesic intersecting S at xe- Since and ur are real near E, the eigenvalue 
splitting is given by the usual formula 

(3) E+-E-= Ah" I iii(0, a;2)5^(0, X2)dx2 + 0(e-(^°+^«)/") 

where cr C S is a neighborhood oi xe- Denote by — 5*^ the phase that comes 
up in (3), where and S*^ stand for suitable Fy{x) as above. By the remark 
following Proposition 1, 5'^ — S*j^ has a non degenerate critical point precisely 
at Xe — n (J. Moreover, the asymptotics of the quasi- modes near Ul{E) 
given in [KaRo], propagate all along T^;, so the integral can be computed by 
standard stationary phase expansion around x = xe- Since the amplitude of 
UR (and ul) is non vanishing, _E+ — E^ is exactly of the order e~So{E)/h^ 

b) The general case. 

Given {pl,Pr) in correspondance, we want to compare, for x € S, the 
action Fy{x) along the bicharacteristic connecting p^ and pR with Agmon 
distance dE'{x) relative to a nearby energy value £". For this, let Te' be a 
minimal geodesic between Ul{E') and Ur{E'), intersecting E at Xe', the left 
and right components of OUe' at = Vl{E') and = yR{E') respectively, 
and consider the lattice of umhilics carrying quasi-modes around such a point. 
At first approximation, umbilics are of the form y = (A^^V2Aiti, A^^-\/2A2t2), 
or by (1), ?/ = {X[^\/2hXiai, X^"^ ^/ 2h\2a2) ^ so the typical neighboring distance 
between umbilics is h{{aih)^^ + (02/1)^^)^^^, which is greater than h/ p, but of 
the same order when tori are not "too flat" . 
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Let y be such an umbilic, and = (Af')^ the corresponding Lagrangian 
manifold. It is easy to see that there is a bicharacteristic 7 C A/, n (and 
points in correspondance) such that Sl — has a non degenerate critical point 
^ e 7 n S. We have : 

{SL-S*ji)^So{E')^2{F^{x)-dE{x))+2{dE{x)-dE'{x))+2{dE'{^^ 

and combining Propositions 1-3 gives, under the above hypotheses, that {Sl — 
5^) - S'o(S') = 0(1), as /i ^ 0, either in case of Birkhoff or KAM tori. To 
compute (3), one has also to know something about the amplitude, so we need to 
improve somewhat the expansions of [KaRo] when getting closer to the caustics. 

References 

[Cr] S.C. Creagh. Tunneling in two dimensions. Proc. on "Tunneling in Com- 
plex systems" (INT 97-1) Seattle, April 30-May 30, 1997. 

[DoSh] S. Dobrokhotov, A. Shafarevich. Math. Phys. Anal. Geom. (2) p.l41- 
177, 1999. 

[Gr] A. Grigis. Seminaire EDP, Expose XXIII. Ecole Poly technique. 1994-95. 
[HeSj] B. Helffer, J. Sjostrand. Comm. Part. Diff. Eq. 9(4) p.337-408, 1984. 
[KaRo] N. Kaidi, M. Rouleux. Comm. Part. Diff. Eq., 27(9 and 10), p.l695- 
1750, 2002. 

[Ma] A. Martinez. Bull. Soc. Math. France 116 (2), p.199-219, 1988. 

[Po] G. Popov. Ann. H. Poincare, Phys. Th. 1(2), p.223-248, 2000. Ann. H. 

Poincare, Phys. Th. 1(2), p.249-279, 2000. 

[Sj] J. Sjostrand. Asymptotic Analysis, 6, p. 29-43, 1992. 

[Wi] M. Wilkinson. Physica 21D, p.341-354, 1986. 



6 



